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Theory of single-molecule experiments in the overstretching force regime 
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We present a statistical mechanics analysis of the finite-size elasticity of biopolymers, consist- 
ing of domains which can exhibit transitions between more than one stable state at large applied 
force. The constant-force (Gibbs) and constant-displacement (Helmholtz) formulations of single 
molecule stretching experiments are shown to converge in the thermodynamic limit. Monte Carlo 
simulations of continuous three dimensional polymers of variable length are carried out, based on 
this formulation. We demonstrate that the experimental force-extension curves for short and long 
chain polymers are described by a unique universal model, despite the differences in chemistry and 
rate-dependence of transition forces. 

PACS numbers: 87.15.By, 83.10.Nn, 87.15.He 



charide such as dextran 
a DNA or RNA strand 



Dynamic force spectroscopy by means of the atomic- 
force microscope (AFM) , laser- or magnetic- tweezers ap- 
paratus, or the biomembrane force probe, allows the di- 
rect probing of the elasticity of individual molecules, and 
as such has rapidly become a mainstay of biophysical re- 
search . These mechanical devices are quite different 
from one another, one prominent difference being their 
equivalent stiffness, in the range of 10~ 4 — lpN/nm for 
tweezers, vs. 10 - 10 5 pN/nm for the AFM [J. The typ- 
ical experiment is a mechanically-induced unfolding of a 
biological polymer made of N domains, e.g. a polysac- 

a protein such as titin Q, 
and so on. As a function 
of increasing force levels different mechanical response 
regimes are observed, beginning with the entropic unfold- 
ing of the polymer chain (now well understood in terms 
of simple worm-like chain (WLC) or freely-jointed chain 
(FJC) models @); to the linear-elastic extension of the 
straightened chain; to the so-called overstretching, typ- 
ically interpreted as a conformational transformation of 
the domain geometry; up to the eventual fracturing of 
the polymer. 

In this Letter we provide a robust statistical mechan- 
ics foundation to the interpretation of the overstretching 
regime, which we describe in terms of the internal dy- 
namics of a chain of two-state systems undergoing a con- 
formational transformation, as described by the double- 
well potential in Fig[TJ For the sake of argument we 
call 'folded' and 'unfolded' the two conformations; how- 
ever the transformation occurs, more generally, between 
two principal local minima of the domain free-energy hy- 
persurface (e.g., for DNA it could as well represent the 
melting transition [ill ] ) . Then, we firstly develop a theo- 
retical model describing experiments at constant applied 
force (a realization of Gibbs ensemble statistics) and we 
show that the conformational change must occur simul- 
taneously for all the domains at a given threshold force. 

On the other hand, experiments performed at 



constant-displacement are a realization of the Helmholtz 
ensemble of statistical mechanics. In our previous work 
Q , we showed that the outcome of the two types of ex- 
periment converge in the thermodynamic limit of infinite 
chain length, N —> oo. In practice, real experiments al- 
ways fall inbetween these two ideal extremes. Therefore, 
here we focus on the intermediate cases described by fi- 
nite values of the k c /k ratio, k and k c being the equivalent 
spring constant (i.e., stiffness) of the domain and of the 
pulling device, respectively. We demonstrate by means 
of Monte Carlo simulations that the typical 'sawtooth' 
pattern Q, observed for the unfolding of large protein 
domains (such as the Ig units in titin), and the 'plateau' 
or kink [H, 0] , observed in the overstretching of DNA and 
polysaccharides (e.g. dextran), have a common origin 
in the size-dependence of the polymer response to the 
external force, the plateau shape being attained in the 
limit of large N. On the same grounds, at a fixed num- 




FIG. 1: Potential energy function with an energy bar- 
rier. Folded and unfolded configurations of the domains are 
schematically represented. 
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ber N of domains, the transition from the 'plateau' to 
the 'sawtooth' response is recovered for increasing values 
of k c /k. Notably, such a behavior of the force-extension 
curves is universal with respect to the specification of 
any additional parameters, such as chemical, structural 
or mechanical constants of the domains. 

We work out a simple model containing the minimal 
ingredients fully describing the overall complex behavior 
of a polymer chain. It consists of an TV-domain, non- 
branched chain clamped at one end, able to describe con- 
formational transitions across an energy barrier. The in- 
ternal state of each domain is described by a potential 
energy V(x) which exhibits two minima corresponding 
to the lengths x = Xf (folded conformation) and x = x u 
(unfolded conformation), connected via an energy barrier 
M at x = xo (see Fig[T|). The energy is written as a C 2 
piecewise function, constructed by imposing continuity 
and differentiability at the joining points x\ and x 2 : 

{\k(x — Xf) 2 < x < xi 

-\k(x - x ) 2 + M X!<x<x 2 (1) 
\k{x — x u ) 2 + AE x>x 2 

For chosen values of the lengths x / and x u , the domain 
spring constant k, and the energy difference AE between 
the two conformations, the other parameters are simply 
given by: S (x u + x f )/2 + 2AE/(kS), 

M = (k/A)[8/2 + 2AE/{k5)] 2 , x x = x f + 6/4 + AE/(kS) 
and X2 = x u — 5/ A + AE/(kS). Therefore, this model 
properly gives a barrier with Xf < xq < x u only for 
\AE\ < kS 2 /4. 

Upon application of a constant force / to the end 
of the polymer identified by the position vector fjv = 
(xat, 1/n, zn) (the other end being fixed in the origin), 
the statistics of the fluctuating chain is a realization of 
the Gibbs ensemble 0. The ensemble partition func- 
tion is given by Z f (f,T) = J J Tn e^l^dq** dp N , with 

Tn = 5R 6W . The augmented Hamiltonian h includes the 
classical kinetic energy of the domains with mass to, their 
total potential energy, and a term, —fzjy, describing the 
applied force along the z-axis Q. In the framework of 
the present minimal model, the partition function can be 
explicitly calculated as: 



Z ' (/,T)= \~) {») W 

x {n (J3k, Pf,x f ,0, x x ) + e' m U (-pk, 13 f, x , x u x 2 ) 
+e- f3AE U((3k,/3f,x u ,x 2 ,+^)} N 
with /3=(fc B T)- 1 and 

f b 2 
II (a, 7, xo, a, 6) = 2 / xe'^^'^ sinh (jx)dx (3) 




The extension r at a given force is obtained from the 
partition function as r = ksT(d\ogZf/df). Since the 
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FIG. 2: (color online) Force-extension curves for the Gibbs 
ensemble: normalized force f/fp vs normalized extension 
r/(Nxf). The black solid lines correspond to different values 
of the energy AE=0, 10, 20, 30, 40, 50 k B T (increasing val- 
ues from the bottom up) for a fixed spring constant k — 2000 
fcflT/(nm) 2 . The blue dashed lines correspond to different 
values of the spring constant fc=10, 15, 30, 100 fcsT/(nm) 2 
(increasing values from the right to the left) for a fixed value 
of the energy barrier AE = 30 fcsT. 



extension is linearly dependent on N, the data for chains 
of different lengths can be scaled to a single curve upon 
diving by N. 

Figure [5] shows the results of the normalized force- 
extension curves, f/fp (where fp 1 = fix;) in terms of 
r/(Nxf) for different values of the energy barrier AE=0, 
10, 20, 30, 40, 50 k B T (black solid lines) at a fixed value 
of fc=2000 fcsT/(nm) 2 , and for different values of the 
spring constant k (blue dashed lines) at a fixed value of 
A_E=30 k^T . Both sets of curves display a force plateau 
at / ~ AE/6, for any AE > 0, with a normalized width 
equal to 6. In our model, the plateau indicates a tran- 
sition in the polymer conformation, meaning that for 
/ < AE/6 each domain is found in the folded confor- 
mation at x=Xf, while for / > AE/S domains are in 
the unfolded conformation at e., the ensemble 

of domains respond cooperatively to the external force. 
Notably, the value of the plateau force inducing the con- 
formation transition does not depend on the spring con- 
stant, k, nor on the temperature. Such a result is readily 
interpreted in the framework of the Bell-Kramers theory 
[5~oj | , as the threshold value of force necessary to make the 
unfolding rate equal to the (reverse) folding one, i.e. low- 
ering the difference AE to zero. For example, in the case 
of dsDNA, which displays a plateau at /=65 pN with a 
6 w 2A, our criterion gives a value AE=3.5 knT, in fair 
agreement with available experimental data [HI HH . A 
similar plateau was observed for other long chain poly- 
mers, such as dextran with N=275, Xf=0.5 nm, x u =0.56 
nm, AE = 13.2 k B T [|, for which the simple criterion 
/ w AE/S gives plateau forces in the range of «900 pN, 
as indeed observed 

While Gibbs ensemble statistics are sampled with a 
constant applied force, a dual situation can be realized 
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FIG. 3: (color online) Monte Carlo force-extension curves at 
T—293 K, for: (top panel) different decreasing values of the 
device spring constant k c =5, 2, 1, 0.5, 0.01 fcsT/ (nm) 2 (from 
the top down) and N=i; (lower panel) increasing number of 
domains N = 4,30,300 with k c = 2 fc B T/(nm) 2 (bottom 
panel). The red dashed line corresponds to the Gibbs ensem- 
ble. The remaining parameters are AE = 30 fcsT, Xf = 2.5 
nm, x u — 3xf and k = 100 feT/nm 2 . 



by imposing the extension, i.e. by controlling the poly- 
mer end-to-end distance. The statistics of the fluctu- 
ating polymer in this latter scheme is a realization of 
the Helmholtz ensemble. As shown in Ref.0], the cor- 
responding partition function Z r cannot be written in 
closed form and, as opposed to the Gibbs case, the cor- 
responding extension r is non-linearly dependent on N. 
However, we showed that the partition functions in the 
two ensembles are formally related via a Laplace trans- 
form, and we demonstrated [|[ that they lead to a com- 
mon force-extension curve in the thermodynamic limit. 

It should be noted that any AFM or tweezers exper- 
iment falls in an intermediate regime between the two 
ideal extremes, of purely constant-force or constant- 
extension, since either constraint on the terminal domain 
of the chain is mediated by a mechanical device (such 
as the AFM cantilever, or the laser-bound microsphere, 
plus a molecular spacer providing adhesion). The device 
is characterized by its own effective elastic constant k c , 
which is coupled in series to the chain of domain springs 
k. In the limit of a soft device, k c /k — >0, the statistics of 
the coupled system reduces to the Gibbs ensemble for the 
isolated molecule fluctuating under a constant force. On 



the other hand, for a very stiff device, k c /k — > oo, one re- 
covers the Helmholtz ensemble for the isolated molecule 
held at a fixed extension by the fluctuating force (l4j |. 
To describe such a situation, we adopt a Monte Carlo 
(MC) numerical approach, simulating the stretching of 
the chain produced by a device with a proper adjustable 
elastic stiffness. 

In Figure El top panel, we report the results of the MC 
simulations at T=293 K, for decreasing values of the k c /k 
ratio, from 0.05, that is well within the Helmholtz statis- 
tics regime, down to 1 x 10~ 4 , i.e., approaching Gibbs 
ensemble statistics. The remaining parameters are set 
to N=A, AE=30 k B T, xf=2.5 nm, x u =3x f and fc=100 
fcsT/nm 2 , which can be considered representative of a 
medium-sized, multi-domain chain protein. At large val- 
ues of k c /k, the domains exhibit a sequence of indepen- 
dent conformational transitions to the unfolded configu- 
ration, generating a series of TV peaks (sawtooth pattern) 
which closely resemble the experimental results obtained 
for short chains (e.g., a titin fragment with 7V=8, Xf=4 
nm, x u =32 nm, A£=ll.l k B T [l|). For k c /k -> 
the peak-to-valley width, A/, of the sawtooth shrinks 
and the curve approaches the k c =0 cooperative plateau 
of Gibbs statistics. In substantial agreement with this 
finding, pulling experiments on native titin by means of 
optical tweezers [15| . having a very small equivalent k c 
compared to the AFM one, do not reveal the sawtooth 
pattern, but rather a smooth, monotonic branch reminis- 
cent of the horizontal plateau. 

On the other hand, a similar asymptotic trend is ob- 
served (Fig. [3j bottom panel) when the chain length, i.e. 
the number of domains, is increased, at a fixed value of 
k c /k. As N increases, the width A/ is decreased un- 
til, at a large enough N, the force-extension curves ap- 
proach again the plateau curve of the Gibbs ensemble. 
It is worth noting that a similar trend was observed in 
experiments performed on native titin, comprising sev- 
eral hundreds of Ig domains, for which the width A/ was 
of the order of 80 pN Q , compared to the much shorter 
8-monomer titin, for which A/ > 200 pN. The experi- 
ments performed on dextran, a long polysaccharide with 
N = 275, [H, HH whose response to the applied force 
shows a plateau closer to the typical DNA-like behav- 
ior, can also be rationalized on this basis. In summary, 
we proved that the macroscopically different behavior 
of small-TV polymers (such as titin) vs. long polymers 
(such as dextran, DNA), as well as experiments done on 
a same polymer but with devices having widely different 
stiffness, can be interpreted with the very same unifying 
model, interpolating between the two extremes of pure 
Gibbs or Helmholtz statistics. 

As observed by several authors, each branch of the 
sawtooth pattern can be nicely fitted by a sequence of 
FJC, or WLC curves (see FigJH left panel, dashed lines) 
with a proper value of the persistence length, up to the 
unfolding of each domain (see e.g. titin spectrin flij ]. 
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FIG. 4: (color online) Left panel: definition of A/ and 6 for a typical force-extension curve with N = 4. Dashed lines for 
the growing branches fitted to FJC with increasing contour length. Central panel: plot of log(A///^) vs. log(JV) for k c —l 
fcsT/(nm) 2 and k/k c =20, 30, 100 (blue, red and black line, respectively). Right panel: plot of tan(#) vs. k c for iV=4, and the 
same k values. Remaining parameters AE — 30 fcsT, xj = 2.5 nm, x u — 3xf and T—293 K. 



fibronectin |17| . synaptotagmin 18 1). Beyond this point, 
the force relaxes to a smaller value, until the next curve is 
met and the force can start rising again upon increasing 
displacement. 

Since the physical origin of the growing branch of the 
curves is well understood on the basis of FJC or WLC 
models, we analyzed the decreasing branch, as identi- 
fied by the common width A/ and angle 8 in Fig0] (left 
panel) which were extracted from our MC simulations as 
a function of TV and k c /k. 

By looking at Fig. [4] (center), the peak-to- valley 
width shows a power-law decrease with the chain length, 
A/ ~ N~ a , the exponent a — 1.3 being remarkably 
independent on the k c /k ratio. This finding indicates 
that attainment of the thermodynamic limit is mainly 
dictated by the thermal force scale, fp, and to a much 
lesser extent by other structural and chemical details of 
the polymer. It is worth noting that the value of the 
exponent is in agreement with previous results on mono- 
stable FJC and WLC models with extensible bonds 

The last plot on the right of Fig. @] reports the behav- 
ior of tan(0) as a function of the device stiffness, k c . The 
observed linear dependence is another remarkable result, 
completely describing the transition between the two ex- 
tremes (Gibbs and Hclmholtz ensembles), while taking 
into account all the intermediate cases. For k c /k — > oo 
we have tan(0) — > — oo or, cquivalently, 8 n/2. 
In other words, the decreasing branches of the force- 
extension curve must be exactly vertical in the case of 
the Hclmholtz ensemble. Notably, both the values of A/ 
and 8 are fully prescribed, i.e. the entire shape of the 
force-extension curve is uniquely defined, once the free 
parameters of the model are specified. 

In conclusion, wc described the statistical mechanics 
of chain polymers composed by domains with two stable 
states, subject to a pulling force by a molecular-scale me- 
chanical device. We showed that for short chain length, 
or large stiffness of the device, the domain response is 



uncorrelated and originates the typical sawtooth force- 
extension curve observed in many experiments. On the 
other hand, upon increasing chain length, or vanishing 
device stiffness, the response is cooperative and results 
in the plateau-like curve, also observed in other experi- 
ments. Despite the simplicity of the model, such a frame- 
work provides a unified picture for such apparently con- 
trasting experimental situations. 

F.M. acknowledges the University of Cagliari for the 
extended visiting grant and the IEMN, University of Lille 
I, for the kind hospitality. L.C. acknowledges financial 
support by RAS under the project 'Modellizzazionc Mul- 
tiscala della Meccanica dei Materiali Complessi (M4C)'. 



[1 
[2 
[3 
[4 
[5 
[6 
[7 

[8 

[9 
[10 

[11 
[12 

[13 

[14 

[15 
[16 

fi- 
lls 



Ritort, J. Phys.: Condens. Matter 18, R531 (2006) 
Cleri, Sci. Model. Simul. 15, 369 (2008) 
R. Chaurasiya et al. Phys. Life Rev. 7, 299 (2010) 
C. Neuman, A. Nagy, Nature Meth. 5, 491 (2008) 
Rief et al, Science 275, 28 (1997). 

Science 276 1109 (1997). 
Y. Cui, C. Bustamante, Science 271 795 



E. D. Siggia, Macromolecules 28, 8759 



Rief et al.. 
M. Smith, 
(1996) 

J. F. Marko, 
(1995). 

F. Manca et al., J. Chem. Phys. 136, 154906 (2012) 

G. I. Bell, Science 200, 618 (1978) 

P. Cluzel et al., Science 271, 792 (1997) 

A. Ahsan, J. Rudnick, R. Bruinsma, Biophys. J. 74, 132 

(1998) 

M. Rief, J. M. Fernandez, H. E. Gaub, Phys. Rev. Lett. 
81 4764 (1998) 

H. J. Kreuzer, S. H. Payne, Phys. Rev. E 63, 021906 
(2001) 

M. S. Kellermayer et al., Science 276, 1112 (1997) 
M. Rief et al., J. Mol. Biol. 289 553 (1999) 
A. F. Oberhauser et al., J. Mol. Biol. 319 433 (2002) 
M. Carrion- Vazquez et al., Prog. Biophys. Mol. Biol. 74 
63 (2000) 



